Surface contribution to the anisotropy of magnetic nanoparticles 
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We calculate the contribution of the Neel surface anisotropy to the effective anisotropy of magnetic 
nanoparticles of spherical shape cut out of a simple cubic lattice. The effective anisotropy arises 
because deviations of atomic magnetizations from collinearity and thus the energy depends on the 
orientation of the global magnetization. The result is second order in the Neel surface anisotropy, 
scales with the particle's volume and has cubic symmetry with preferred directions [±f , ±1, ±f]. 
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With the decreasing size of magnetic particles, sur- 
face effects are believed to become more and more pro- 
nounced. A simple argument based on the estimation 
of the fraction of surface atoms shows that for a parti- 
cle of spherical shape and diameter D (in units of the 
lattice spacing), this fraction is an appreciable number 
of order 6/D. Regarding the fundamental property of 
magnetic particles, the magnetic anisotropy, the role of 
surface atoms is augmented by the fact that these atoms 
in many cases experience surface anisotropy (SA) that 
by far exceeds the bulk anisotropy. As was suggested by 
Neel |j| and microscopically shown in Ref . , the leading 
contribution to the anisotropy is due to pairs of atoms 
and can be written as 
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where rrij is the reduced magnetization (spin polariza- 
tion) of the ith atom, are unit vectors directed from 
the itli atom to its neighbors, and Ly is the pair- 
anisotropy coupling that depends on the distance be- 
tween atoms. Eq. ([!]) describes in a unique form both the 
bulk anisotropy including the effect of elastic strains and 
the effect of missing neighbors at the surface that leads to 
the SA. In particular, for an unstrained simple cubic (sc) 
lattice the bulk anisotropy in Eq. ([!]) disappears since 
m 2 + iriy + mr z = 1 is an irrelevant constant, and one has 
to take into account the dropped (much smaller) terms 
of Eq. (jl]) that yield the cubic bulk anisotropy. On the 
other hand, surface atoms experience (large) anisotropy 
of order L due to the broken symmetry of their crys- 
tal environment - the so-called Neel surface anisotropy 
(NSA) . These atoms can make a contribution to the effec- 
tive volume anisotropy decreasing as 1/D with the par- 
ticle's linear size: Kv, e s = Kv + Ks/D, as was observed 
in a number of experiments (see, e.g., Refs. || Q). 

The 1/D surface contribution to Ky. c s is in accord 
with the picture of all magnetic atoms tightly bound 
by the exchange interaction whereas only the surface 
atoms feel the surface anisotropy. This is definitely true 
for magnetic films where a huge surface contribution to 
the effective anisotropy has been observed. The same 
is the case for cobalt nanoclusters of the form of trun- 
cated octahedrons H where contributions from different 
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FIG. 1: Magnetic structure of a spherical nanoparticle of 
linear size N — 15 with L/J = 2 for the global magnetization 
directed along [1,1,0], showing atoms in the plane 2 = 0. 



faces, edges, and apexes compete resulting in a nonzero, 
although significantly reduced, surface contribution to 
Kv, e {f- However, for symmetric particle shapes such as 
cubes or spheres, the symmetry leads to vanishing of this 
(first-order) contribution. In this case one has to take 
into account deviations from the collinearity of atomic 
spins that result from the competition of the SA and the 
exchange interaction J. The resulting structures (for the 
simplified radial SA model) can be found in Refs. |6|, ^, ||] 
(see also Fig. [l] for the NSA). In the case L > J devia- 
tions from collinearity are very strong, and it is difficult 
if not impossible to characterize the particle by a global 
magnetization suitable for the definition of the effective 
anisotropy. On the other hand, in the typical case L< J 
the magnetic structure is nearly collinear with small devi- 
ations that can be computed perturbatively in L/J <C 1. 
The global magnetization vector mo can be used to de- 
fine the anisotropic energy of the whole particle. The key 
point is that deviations from collinearity and thus the en- 
ergies of the system are different for different orientations 
of mo, even for a particle of a spherical shape, due to the 
crystal lattice. For the latter the overall anisotropy per 
unit cell is proportional to L 2 /J, i.e., it scales with the 
particle's volume. 



The aim of this Letter is to illustrate this idea by cal- 
culating the second-order contribution from the NSA to 
the effective particle's anisotropy for the minimal model 
of a magnetic nanoparticle of spherical shape cut out of a 
sc lattice. We will neglect the small cubic anisotropy and 
magnetostatic effects for the sake of transparency. The 
problem will be solved numerically on the lattice by min- 
imizing the energy with the help of a damped Landau- 
Lifshitz equation without the precession term, with the 
average particle's magnetization constrained in a desired 
direction. We also produce an analytical solution in the 
continuous limit of larger particles that will be shown to 
agree with the numerical solution. 

We consider the nearest-neighbor form of Eq. (Q) with 
the unique constant L. For a sc lattice it reduces to 
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where Zi a = 0,1,2 are the numbers of available near- 
est neighbors of the atom i along the axis a. One can 
see that the NSA is in general biaxial. For L > and 
Zi a — < Zip = 1 < Zi-y — 2 the a-axis is the easy axis 
and the 7-axis is the hard axis. If the local magnetiza- 
tions m.i are all directed along one of the crystallographic 
axes a, then the anisotropy fields Hai — —dHAi/dnii 
are also directed along a and are thus collinear with m.;. 
Hence, at least for L « J, there are no deviations from 
collinearity if the global magnetization mo is directed 
along one of the crystallographic axes. For other orien- 
tations of mo, the vectors m, and H^j are not collinear, 
and the transverse component of Hai with respect to 
mj causes a slight canting of m 4 and thereby a devia- 
tion from the collinearity of magnetizations on different 
sites. This adjustment of the magnetization to the sur- 
face anisotropy leads to the lowering of energy. As we 
shall see, this effect is strongest for the [±1,±1,±1] ori- 
entations of mo. For both signs of L these are easy ori- 
entations, whereas [±1,0,0], [0, ±1,0], and [0,0, ±1] are 
hard orientations. 

We consider here explicitly spherical particles cut out 
of a cube with dimensions N x N x N in the units of 
the atomic spacing. If an atom is within or exactly on 
the sphere with the diameter D = N — 1, it belongs 
to the particle. The number of atoms in the particle 
J\f approaches M ^ (tt/6)(./V - l) 3 for N > 10, with 
fluctuations for smaller N. Our numerical results for the 
magnetic energy of spherical particles as a function of 
the orientation of the global (average) magnetization are 
shown in Fig. |[ They confirm the statements of the 
previous paragraph. 

To produce Fig. 0, we use the classical Hamiltonian 
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with the nearest-neighbor exchange coupling J and Ha of 
Eq. (0). To fix the global magnetization of the particle 




FIG. 2: Reduced shifted energy of the particle for different 
orientations of its global magnetization, obtained from Eq. 
(^) . These curves manifest the cubic symmetry of the effective 
anisotropy, see Eq. ([jl]). 



in a desired direction Vq (uq\ = 1), we use the energy 
function with a Lagrange multiplier A: 
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To minimize we solve the evolution equations 



rhi = - [mi x [mi x Fj]] , 
A = dF/d\=-N'(v-v ) 



-dT/dnii 
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starting from mi = Uq = m and A = 0, until a station- 
ary state is reached. In this state v = v$ and [m* x Fj] = 
0, i.e., the torque due to the term J\[\-(v — Vq) in T 
compensates for the torque acting to rotate the global 
magnetizations towards the minimum-energy directions 
[±1,±1,±1]. Since the former torque is unphysical, this 
method is applicable only for a small surface anisotropy, 
so that both torques are small, and adding a small for- 
mal compensative torque does not strongly distort the 
magnetic structure. 

In physical terms, the existence of the well-defined 
state with a given orientation of the global magnetiza- 
tion can be justified as follows. For L <C J, the relax- 
ation of the magnetization splits into two stages. The 
first stage, adjustment of the magnetic structure to the 
surface anisotropy, involves energies of order L and is 
relatively fast. The second stage, rotation of the global 
magnetization to the global energy minimum with the 
magnetic structure adjusted at any moment, involves en- 
ergies of order L 2 /J and is much slower. Introducing the 
global-orientation constraint above eliminates the second 
stage of the relaxation, so that the result of the first re- 
laxation stage is seen in pure form. 

Fig. H shows the dependence of the normalized parti- 
cle energy differences between the basic directions [001], 
[011], and [111]. One can see that AE/Af tends to a 
large-A limit, i.e., for large linear sizes N the energy 
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FIG. 3: Differences of the particle energies between main 
orientations of the global magnetization vs the particle size 
in the scaled form for L/J = 0.1 and 0.01. The scaling is 
valid for TV < J/L, and its violation for L/J = 0.1 is seen in 
the right part of the figure. 



differences due to the SA scale with particle's volume 
V oc W ~ (N - l) 3 . These results suggest that the 
problem can be solved analytically with the help of the 
continuous approximation for N 3> 1. To this end, we 
replace in Eq. (^) the number of nearest neighbors of a 
surface atom by its average value 

Zia z ia =2 - \n a \/ max {\n x \, \n v \, \n z \} . (6) 

Here n a is the a-component of the normal to the surface 
n. The surface-energy density can then be obtained by 
dropping the constant term and multiplying Eq. (^) by 
the surface atomic density /(n) = max{|naj|, \n y \, \n z \}: 

E s (m,n) = -- [\n x \m 2 x + \n y \m 2 y + \n z \m 2 z ] . (7) 

At equilibrium, in the continuous approximation the 
Landau-Lifshtz equation reduces to 

m x H off = 0, H off = + JAm, (8) 

where A is the Laplace operator and the anisotropy field 
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For L <C J the deviations of m(r) from the homogeneous 
state mo are small and one can linearize the problem: 



m(r) = m + ip(r, m ), ip = \ip\ <C 1. 



(10) 



The correction ip is the solution of the internal Neumann 
boundary problem for a sphere 



or 
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where n = r/R and m stands for mo with the index 
dropped for transparency, ip has the form 

tp(r,m) = — / d 2 r'G(r,r')f(m,n') (12) 
47r Js 



with the Green function 
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S(r,r') = ^R 4 + r 2 r' 2 ~2R 2 (r-r>). (13) 

One can make the estimation 

ip(r, m) ~ RL/J ~ NL/J. |r| = R. (14) 

This shows that for whatever small values of L the appli- 
cability condition of our linearization method ip = \ip\ <C 
1 will be invalidated for sufficiently large particle sizes. 

Now we are prepared to calculate the magnetic energy 
of the nanoparticle. Dropping the trivial constant term 
leads to the second-order energy 

fta .^ + ^_jr (( ^ W )' + /A(^.*) 

(15) 

that is a sum of the inhomogeneous exchange and 
anisotropy energies. With the help of Eq. ( pM| ) this yields 

E2 ~hh\ Is d2rdVG ( r ' r >K n ' n ') ( 16 ) 
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The first term in <i>(m, n, n') can be simplified using 
m • dEs(m, n)/dm = 2Es(m, n) following from Eq. (Q). 
The second term in $(m, n, n') is quadratic in the mag- 
netization components and contributes only with the ir- 
relevant term proportional to m 2 + m 2 + m 2 = 1 to the 
energy. Thus £2 simplifies to 
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that is of fourth order in the global-magnetization com- 
ponents m a . Taking into account the cubic symmetry 
and computing numerically a double surface integral one 
can write the result of Eq. ( |l8| ) as 
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where Jo — zJ = 6J. This defines the large-iV asymp- 
totes in Fig. H that are shown by the horizontal lines. 

The analytical results above are valid for particle sizes 
N in the range 



1< N « J/L. 



(20) 
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The lower boundary is the applicability condition of 
the continuous approximation. Since the surface of a 
nanoparticle is made of atomic terraces separated by 
atomic steps, each terrace and each step with its own 
form of NS A [see Eq. (||)] , the variation of the local NSA 
along the surface is very strong. Approximating this vari- 
ation by a continuous function according to Eq. ([}]) re- 
quires pretty large particle sizes N. This is manifested by 
a slow convergence to the large- N results in Fig. ||. 

The upper boundary in Eq. ( pQ ) is the applicability 
condition of the linear approximation in if>, see Eqs. ( |l0| ) 
and (|b|). For N > J/L deviations from the collinear 
state are strong, and the effective anisotropy of a mag- 
netic nanoparticle cannot be introduced. The solution 
found above becomes invalid even for orientations of 
the global magnetization along the crystallographic axes 
where ■0 = 0. In this case those surface spins close to 
the equatorial plane (aim) for L > or to the poles 
(n||m) for L < develop instability and turn away from 
m for N > J/L. Gradual disappearance of the collinear 
magnetic structure of a particle with increasing size stems 
from the "softnening" of the exchange interaction at large 
distances. A related phenomenon is the breakdown of 
the single-domain state of particles with a uniaxial bulk 
anisotropy with increasing size due to the magnetostatic 
effect. 

As we have seen in Eq. ( pj| ) , the contribution of the SA 
into the overall anisotropy of a magnetic particle scales 
with its volume V oc N 3 ~ J\f. This surprising result, 
that contradicts the initial guess on the role of the sur- 
face effects based on the ratio of the numbers of sur- 
face and volume spins ~ 6/D, is due to the penetration 
of perturbations from the surface deeply into the bulk. 
If a uniaxial bulk anisotropy Dy is present in the sys- 
tem, perturbations from the surface will be screened at 
the bulk correlation length (or the domain-wall width) 
8 ~ V 'J/D v . Then for D ~ N > 5 the contribution of 
the SA to the overall anisotropy will scale as the sur- 
face: £ 2 ~ {L 2 /J) N 2 8. As follows from Eq. @, this 
regime requires Dy ^ L 2 / J, i.e., the dominance of the 
bulk anisotropy over the SA in the overall anisotropy. 

In most cases the bulk anisotropy is much smaller than 
the surface anisotropy for the microscopic reasons dis- 
cussed at the beginning of this Letter. Then, at least 
for not too large particles, N < 8, contributions of both 
anisotropics to the overall anisotropy are additive and 
scale as the volume. If the bulk anisotropy is cubic, both 
contributions have the same cubic symmetry [sec Eq. 
(|l9|)1, and the experiment should yield a value of the effec- 
tive cubic anisotropy different form the bulk value For 
the uniaxial bulk anisotropy, the two contributions have 
different functional forms. Even if the bulk anisotropy 
is dominant so that the energy minima are realized for 
m||e z , the surface anisotropy makes the energy depen- 
dent on the azimuthal angle <p. This changes the type of 



the energy barrier for the particle creating saddle points. 
The latter, in particular, strongly influences the process 
of thermal activation of magnetic particles |J. 

We stress that we have calculated the second-order 
contribution of the Neel surface anisotropy to the effec- 
tive anisotropy of a magnetic particle, and this is the 
only effect for symmetric particle shapes such as cubic 
or spherical. For small deviations from this symmetry, 
i.e., for weakly elliptic or weakly rectangular particles, 
there is a correspondingly weak first-order contribution 
£ i that adds up with our second-order contribution. For 
an ellipsoid with axes a and b = a(l + e), e <C 1, one has 
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can be large even for L/J <C 1. Whereas £\ scales with 
the particle's surface and can be experimentally identified 
as a surface contribution, £2 scales with the volume and 
thus renormalizes the volume anisotropy of nanoparticlcs. 

The Neel constant L is in most cases poorly known. 
However, for metallic Co Ref. |l0| quotes the value of 
SA -1.5 x 10 8 erg/cm 3 , i.e., L ~ -10 K. This is much 
smaller than J ~ 10 3 K, which makes our theory valid 
for particle sizes up to N ~ J/L ^100, according to Eq. 
(pp|). For this limiting size one has £n/£\ ~ 1/e that is 
large for nearly spherical particles, e -C 1. 

We are indebted to R. Schilling for critical reading of 
the manuscript. D. G. thanks A. A. Lokshin for a valu- 
able discussion. 
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